Flatband localization in the Anderson-Falicov-Kimball model 
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The insulator-metal-insulator transition caused by a flatband is analyzed within dynamical mean- 
field theory using the Anderson-Falicov-Kimball model. We observe quantitative disagreement be- 
tween the present approach and previous results. The presence of interactions enhances derealiza- 
tion. 
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I. INTRODUCTION 

Recent studies have shown a new disorder-induced 
insulator-metal transition (MIT) of one-electron states, 
called the "inverse Anderson transition" ^ The exis- 
tence of this transition becomes already visible for non- 
interacting particles when the system has highly degen- 
erated localized states forming a flatband. In a flatband 
localized states may melt into extended states due to dis- 
order. For weak disorder, the localization is not a conse- 
quence of the strength of disorder but of the flatband. In- 
creasing the degree of disorder localization-delocalization 
and delocalization-localization transitions appear. 

New aspects must appear when we also consider in- 
teractions between particles of the system. The present 
work addresses this issue. The investigation of this prob- 
lem is particularly rich. The interaction causes a Mott- 
Hubbard metal- insulator transition (MIT) 3 that com- 
petes with the Anderson transition^ 

The inverse Anderson transition has been studied us- 
ing various numerical techniques, among them, level 
statistics, 5 f(a) characteristics of wave functions 6 - and 
the participation radius 7 framework. However, these 
powerful tools are not adequate for interacting systems 
in which matrices of high order must be exactly diago- 
nalized restricting the study to extremely small systems. 

The dynamical mean-field theory (DMFT) 8 is a good 
tool to investigate the Mott-Hubbard-Anderson MIT in 
lattice electrons with local interactions and disorder. The 
Anderson transition has for example been explored on 
the Bethe lattice considering the Hubbard 9 and Falikov- 
Kimball models^ The metal and the insulator phases 
are detected by analyzing directly the local density of 
states (LDOS). The averaged LDOS can vanish in the 
band center at a critical disorder strength for a wide va- 
riety of averagesi 11 ' 12 In particular, the arithmetic mean 
of this random one-particle quantity is non-critical at 
the Anderson transition and hence cannot help to detect 
the localization transition. By contrast, the geometric 
mean gives a better estimate of the averaged value of the 
LDOS, 9,13 as it vanishes at a critical disorder strength 
and hence provides an explicit criterion for the Anderson 
localization ^ 14 : 15 We have adopted the Holder mean and 
analyzed how the averaged LDOS depends on the Holder 



parameter. 

In this paper, we investigate the disorder-induced 
insulator-metal transition in a flatband of the Falikov- 
Kimball model using the DMFT. We present the ground- 
state phase diagram for different values of a parameter 
measuring the disorder strength and the dependence of 
the MIT transition on the Coulomb repulsion between 
fermions. Recently, we have applied the DMFT to this 
model. First, we showed that not only the geometric 
mean can offer a good approximation for the averaged 
LDOS providing an explicit criterion for Anderson local- 
ization. We found that the averaged LDOS can vanish in 
the band center at a critical disorder strength for a wide 
variety of generalized Holder meani^i Second, we have 
analyzed how the presence of the next-nearest-neighbor 
hopping influences the phase diagram of the ground state 
of this model, and, third, we studied the main effects of 
the long-range correlated disorder^ 2 - 

This paper is organized as follows. In the next section 
HTl we introduce the Anderson-Falikov-Kimball model. 
The DMFT approach is described in section IIIII In sec- 
tion IIVI we present the results concerning the phase dia- 
gram. Finally in section [V] we conclude. 



II. MODEL 

The Anderson-Falicov-Kimball modeli 3 - is a tight bind- 
ing model having two species of fermionic particles, mo- 
bile and immobile, which interact with each other when 
both are on the same lattice site. We introduce a local 
random potential for the mobile particles, giving rise to 
a competition between interaction and disorder. This 
model has been applied to mixed valence compounds 
of rare earth, transition metal oxides, binary alloys and 
metal ammonia solutions. 16 Its Hamiltonian is 



(1) 



(>j) 



where cf (cj) and (/*) are, respectively, the cre- 
ation (annihilation) operators for the mobile and immo- 
bile fermions at lattice site i. q is a random potential 
describing the local disorder, t is the electron transfer 



integral connecting nearest-neighbor sites, and U is the 
Coulomb repulsion when mobile and immobile particles 
occupy the same site. We consider that the occupation 
of immobile particles is site independent having a prob- 
ability p = 1/2. The number of mobile particles on site 
i is given by rii = c\ c$. A chemical potential /i is in- 
troduced for the mobile subsystem to fix the system in 
the half-filled band (nj = 1/2). Here, the energy will be 
given in units of the hopping element t (i.e., t = 1). 



III. DYNAMICAL MEAN-FIELD EQUATIONS 

The DMFT is calculated from the Hilbert transform 

du>N (uj) 



G(E) 



(2) 



n(E) -lu + 1/G{E) ' 

where Ao(e) is the non- interacting density of states, G(E) 
the translationally invariant Green function, and r)(E) a 
hybridization function describing the coupling of a se- 
lected lattice site with the rest of the systemJ^ For the 
flatband the non-interacting density of states is 

N (E) = ^[5(E-E ) + 5(E + E )}, (3) 

where the highly degenerated energies are ±Eq>. The re- 
lation between G(E) and rj(E) is obtained in a straight- 
forward way from Eq. ((2|) and Eq. ([3]) 



V (E) = 



i 



v 4G(£) 2 
The LDOS is given by±2 

1 



El- 



2G(E) 



P(E,ei) = Im G(E,ei) , 

7T 



(4) 



(5) 



where G(E,Ci) is the local Q-dependent Green func- 
tion. For the Anderson- Falicov-Kimball model we obtain 



{U/2f 



U 2 o? 



(6) 



7T [a? + s 2 + (C//2) 2 ] 2 

where on — E — e$ — r and r and s are, respectively, the 
real and imaginary parts of T]{E). 

We consider that ej is an independent random variable 
characterized by a probability function p(e^) = $(A/2 — 
|ej|)/A, with $ being the step function. The parameter A 
is a measure for the disorder strength. The self-consistent 
DMFT equations are closed inserting 



where 



G{E) 



duj' 



E 



1/9 



(7) 



(8) 
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FIG. 1: Density of states for disorder strength A = and 
U = 1, 4, 8 and 12. 
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FIG. 2: Averaged local density of states at U — and disorder 
strength A = 1, 5, and 10 for (a) q = 1 and (b) q = 0. 



IV. RESULTS 



The parameter q defines the ^-Holder average. The arith- 
metic and geometric mean are found, respectively, using 
q = 1 and q — > 0. 



First, we investigate the case without disorder. In this 
limit, where A = 0, we find, independently of q, that 
the analytical expression for the density of states of the 
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FIG. 3: Ground-state phase diagram as a function of energy 
for [7 = 0. The circles were taken from Ref. [1]. The triangles 
are determined from the numerical solution of the DMFT 
equations for q = 0. The dashed region is the energy area 
in which we cannot estimate the transition points accurately. 
Lines are guides to the eye. 



FIG. 4: Averaged local density of states at U — 1 and disorder 
strength A = 0.1, 0.5, 1.0 and 5.0 for q = 0. 



flatband is 



P q {E) 



E 



V[# 2 ~ Wo + U/2Y][{E - £7/2)2 - W] 



(9) 

if \{E - U/2)\ < \E\ < \(E + C//2)| and P q (E) = 0, 
otherwise. Typical results are shown in Fig. [T] In order 
to compare our results with Refs. [1] and [2], we have set 
the highly degenerated state at Eq = 4. Below [7 = 8 
we have two bands of bandwidth [7. For [7 > 8 we have 
also two bands, but now the bandwidth is constant and 
equal to 8. The bandgap is always 8 — U\. This bandgap 
for [7 ^> 1 exhibits the Mott insulator relationship be- 
tween the bandgap and the repulsive Coulomb potential 
(bandgap ~ U). Note that for U = 8 the system presents 
only one band of bandwidth 8, and only for this case 
P,(0) + 0. 

Next, let us explore the case A / 0. Here, the results 
are obtained numerically. We considered as initial P q (E) 
a uniform distribution with bandwidth greater than the 
Lifshitz band edge. Then we determined G(E) in order 
to obtain rj(E) and finally the new values of P q (E). This 
procedure is repeated until we find the stable configura- 
tion. 

Fig. [2] shows the energy dependence of the averaged 
LDOS for [7 = and typical values of A. Note that 
the inclusion of the disorder (i.e., A ^ 0) suppresses 
the highly degenerated localized states. We obtain the 
Anderson localization for a fixed q varying the disorder 
strength A for each value of U, and then determining 
the values of A when P q {E) = 0. For the highly de- 
generated energy state, the arithmetic mcan(q = 1) of 
the LDOS does not vanish at a finite critical disorder 
strength. Hence, we consider it to be non-critical at 
the Anderson transition. Using the geometric average 
(q = 0), the LDOS vanishes at the highly degenerated 
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FIG. 5: Averaged local density of states at U — 4 and disorder 
strength A = 1, 2, 5 and 10 for q = 0. 



energy state for a finite value of A. The detection of the 
Anderson localization is obtained using q = 0. 

Fig. [3] presents the phase diagram of the ground state 
for the Anderson-Falikov-Kimball model as a function of 
energy for [7 = using the geometrical average. Our 
results are represented by triangles. The results of Ref. 
[1] are marked by circles. As we use an iterative pro- 
cess, our P q {E) not always converges to a stable value 
for large A. Within the dashed area of Fig. 3 we can- 
not estimate the transition points accurately. We can 
observe three phases: extended gapless, localized gapless 
and band gap.— 

The DMFT gives different results from the ones ob- 
tained in Refs. [1] and [2]. The DMFT approach reduces 
the extended phase, increasing the critical disorder A c for 
the localization-delocalization transition and decreasing 
the one for the delocalization-localization transition. The 
considerable quantitative disagreement between both ap- 
proaches requires new investigations to better understand 
these differences. 
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FIG. 6: Ground-state phase diagram as function of energy 
for U — 0, 1 and 4 obtained from the numerical solution of 
the DMFT equations. The dashed region is the energy area 
in which we cannot estimate the transition points accurately. 
Lines are guides to the eye. 

We now consider the influence of the Coulomb repul- 
sion U on the results. Fig. [4] shows the averaged LDOS 
for A = 0.1, 0.5, 1.0 and 5.0 at U = 1. The results 
for U = 4 are exhibited in Fig. [51 for A = 1, 2, 5 and 
10. As already observed for U = 0, we find two symmet- 
ric bands. However, for small A, the LDOS at the end 
bands, corresponding to the smallest and largest ener- 
gies, are higher than those at the band centers. The band 
centers correspond to the highly degenerated states for 
U = 0. Large A destroy the influence of the flatband and 
favour higher LDOS values at band centers and smaller 
ones at the edge bands. When A is large enough, the 
LDOS vanish and the extended phase disappears. For 



fixed A, the bandwidth grows with increasing U. 

Finally, in Figure[S]we present a complete ground-state 
phase diagram for three different values of U, namely 
U = 0, 1 and 4. With U increases the extended phase, 
showing that the interaction may enhance derealization. 

V. CONCLUSIONS 

In the present paper, we studied the solutions of 
the Anderson-Falicov-Kimball model involving a highly 
degenerated localized states forming a flatband. We 
have shown that the new disorder-induced insulator- 
metal transition of one-electron states, called "inverse 
Anderson transition" can be obtained within dynamical 
mean field theory. However, we observed a considerable 
quantitative disagreement between the present approach 
and the previous results*^ The DMFT reduces the ex- 
tended phase, increasing the critical disorder A c for the 
localization-delocalization transition and decreasing A c 
for the delocalization-localization transition. 

We also studied the dependence of the MIT transition 
on the interaction between particles of the model and 
showed that increasing the interaction parameter reduces 
the extended phase for U = 0. 

As the doping of an impurity in the flatband has the 
same characteristics of the impurity states in the quan- 
tum Hall system^, it would be interesting to do similar 
calculations for such a system. 
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